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Abstract 

We present new mathematical alternatives for explaining rotation curves of spiral 
galaxies in the MOND context. For given total masses, it is shown that various 
mathematical alternatives to MOND, while predicting flat rotation curves for large 
radii (r/r^ 3> 4, where rd is the characteristic radius of the galactic disc), predict 
curves with different peculiar features for smaller radii (0.1 < r/rd C 4)- They are 
thus testable against observational data. 


Key words: Gravitation: phenomenology, galaxies: internal motions 
PACS: 04.90.+e, 95.30.Sf, 98.62.Dm 


1 Introduction 


The first mathematical descriptions of the effects of gravity, made by Galileo 
in his study of the free fall of bodies and by Kepler in his study of planetary 
motions, were purely empirical. Though Newton offered a coherent explanation 
of what was behind the laws governing gravitational effects, it was only with 
Einstein’s General Relativity that we had an apparently complete theory of 
gravity. 

However, at the end of the 20 th century, a new enigma concerning the motion 
of ‘celestial bodies’ emerged, in particular, in studying rotation curves of spiral 
galaxies. While Newton’s law of gravity predicts that the velocity of rotation 
in the interior of a galaxy should fall with increasing distance from the galactic 
center if the observed light traces mass, what is observed is the maintenance 
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of a constant velocity with increasing radius, generating flat rotation curves 

[ 11 - 

Two simple ways of dealing with this problem have been suggested: 

(1) assuming that there is more mass ( i.e ., dark matter) in galaxies than is 
observed; 

(2) modifying the law of gravity. 

While much work has been done in the search for possible particle candidates 
for dark matter [2], very little has been done to explore the possibilities of 
modified gravity laws. Until now, the most popular suggestion for a modified 
gravitational law has been Modified Newtonian Dynamics, or, MOND [3-5]. In 
MOND the acceleration a of a body in an external gravitational held is not ex¬ 
actly equal to the acceleration g^ obtained from the Newtonian gravitational 
force. Mathematically, one can write ag = , where g ( x ) is a dimensionless 

function of the ratio x = a/a 0 of the acceleration a to an empirically deter¬ 
mined constant a 0 . Only in the limit a 3> is Newtonian gravity restored. 
The strongest objection to MOND is that it does not have a relativistic the¬ 
ory supporting it. For recent articles criticizing MOND, see Scott et al. (2001) 
[6] and Aguirre et al. (2001) [7]. For a recent positive review of MOND, see 
Sanders (2001) [8]. 

The objective of this letter is to expand the original MOND proposal by pre¬ 
senting mathematical alternatives for the modified gravitational law. Specifi¬ 
cally, we present several alternative mathematical alternative formulations for 
the dimensionless function g, thus following closer the structure of the pio¬ 
neering work of MOND by Milgrom [3-5]. In the next section we present the 
basics of MOND. Simulated rotation curves for several possible MONDian- 
like functions are given in Section 3. The final section presents some brief 
conclusions and perspectives for future work. 


2 MOND 


As discussed in the introduction, the original MOND proposal uses the relation 


ag = g N 


( 1 ) 


2 



where g n is the usual Newtonian acceleration and fi (x) is a function which 
obeys 


j 1 (x > 1 ) 
| x (i«l) 


( 2 ) 


Therefore, in the limit of large accelerations, x 1, the usual Newtonian 
gravity law is obtained. In the other extreme, i<l, however, we have 

a/j, ~ a 2 /a 0 — — GMr~ 2 . (3) 


Thus, using a = V 2 r 
V ~ (GMa 0 ) 1/4 


1 , where V is the rotation velocity of the galaxy, 


( 4 ) 


which is a constant, as is observed for large galactic radii. 

It is common in the literature ( e.g . [ 6 ], [7]) to use the expression 

1 /9 

n(x) — x (l + x 2 ) . (5) 


This formula, proposed by Milgrom [3-5], has the advantage of being invert¬ 
ible. With it one can solve eq. (1) analytically for the acceleration a and, 
consequently, for the rotation velocity Id as a function of the radius r. How¬ 
ever, other functions are also possible, and are discussed in the next section. 


3 Alternative mathematical formulations of MOND 


In his work on the implications of MOND for galaxies [4], Milgrom used as a 
model for a spiral galaxy of total mass M, a disc of mass M d and a central 
spheroidal bulge of mass M s . The fractional masses for the disc and the spher¬ 
ical bulge are a d = M d /M and a> s = M s /M = 1 — a d , respectively, so that the 
total fractional mass 7 = M (r) /M inside a radius r = sr d is 

7 (s) = a^d 0 ) + (1 - a d ) 7s 0 ) , ( 6 ) 

where [4] 
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( 8 ) 


7s («) 


y/8n 3 b 9 

-7 

m t 


(8.5, bs l ^u 
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and 7 (a, z) = / 0 Z e~H a ~ l dt is the incomplete gamma function, b = 7.66924944 
and mt = 2.4082 x 10 -3 are numerical constants. The dimensionless variable 
s = r/r f i = ur/r s is the ratio of the radius r to the characteristic length ry. 
The ratio of r s to r^, u = r s /rd, is less than unity. The radii r ( ] and r s are 
obtained, in practice, by adjusting the luminosity profiles of the spheroidal and 
disc components, using the empirical law of de Vaucoulers for the spherical 
bulge and an exponential function for the disc. 


Following the MOND proposal, we define 


a — fi' ( y ) g K , 


(9) 


where // is a dimensionless function with a dimensionless argument y = 
r^Ja 0 G~ 1 M~ 1 (r), similar to the /j of Milgrom [3-5] in eq. (5). This new func¬ 
tion \j! (■ y ) is such that 


h (y) 


y (y > i) 
i (y «i) 


( 10 ) 


We investigate the following functions // (y) which obey the constraints of eq. 

( 10 ): 


y'i = Vi + y 2 
lA = y coth y 
A 3 = y (1 - e-y)- 1 
H' 4 = y (coth 3 y - 


( 11 ) 
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The behaviour of each of these functions as a function of y can be seen in the 
expansions [9] 


< 


U i 


/4 


/4 


/4 


= Er= 0 r(n-l/2)(2 v /?n!) 1 (-l)""V n 
= 1 + 2/ 2 /2 - r/78 + 2/716 - 5//7128 + ... 

= l + 22/ 2 Er=i (2 y 2 + k 2 7T 2 )~ 1 
~ 1 + 2/ 2 /3 — 2/ 4 /45 + 22/ 6 /245 + ... ( 2 / < n) 

= y E? =0 e~ ky 

~ 1 + y/2 + 2/712 - 2/7720 + 2/730240 + ... (y < 2vr) 
= g[E£i (92/ 2 + ^TT 2 )- 1 ]" 1 
~ 1 + 272/ 2 /45 - (272//45) 4 + (272//45) 6 + ... ( 2 / « 1) 


( 12 ) 


The functions are plotted in Figure 1. 

Using these functions, together with equations ( 6 ), (7) and ( 8 ), we obtain 
curves for the dimensionless rotation velocity v = V (r) / (GMao ) 1 ^ 4 as a func¬ 
tion of s = r/r d for different values of M, a d = M d /M, and u = r s /r d . The 
curves are shown in Figures 2 and 3. 


4 Conclusion 


Inspection of Figures 2 and 3 shows clearly that all the functions y! ( y ) produce 
flat rotation curves. This is true not only for the particular values of M, a d , 
and u of the figures, but for the entire range of physically reasonable values 
for these parameters. Figure 3 shows that a comparison between the curves 
obtained, using the different y! functions presented, together with the original 
Milgrom proposal (eq. (5)), may be useful to distinguish between them, since 
each curve has a peculiar feature in the region 0.1 < r/r d < 4. 

It would be interesting to test the formulas presented here against observa¬ 
tional data, noting that ad and u are not free parameters, but are given by 
the luminosity profiles of the galaxies. The mass M and the constant a 0 are 
the only free parameters to be adjusted. The study of different galaxies gives 
a single value for a 0 , the mass M and the mass-luminosity ratio, M/L, of each 
galaxy, y and y! can lead to different relativistic extensions of MOND, impor¬ 
tant for future studies. For instance, using the expression for the gravitational 
potential ip (r) = — / a ( r) dr, valid for purely radial forces, one can naively 
ascribe a y (r) to the modified gravitational laws obtained with y \, y’ 2 and y' 3 , 
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for example, 


<fVi ( r ) = ““tV 1 + + V GMa o arcsinh ( r \/^) + ’ ( 13 ) 

A /4 (r) = - C ^- + 2 yjGMa 0 f] — arctan 7 ) + Ao , (14) 

2 r v ^ ri7r \rnr\GMJ 


and 


^ (r) = [ ln ( r Vill) + £ K {~ nr \fdK) 


+ Ao j 


(15) 


where (40 is a constant of integration. 

Therefore, in the search for a complete theory for MOND, it is important to 
study alternative MONDian functions. The MONDian functions given in this 
letter can be seen as a step in this direction. 
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Fig. 1. Top: Curves showing the behaviour of the y! functions as a function 
y = ry / aoG l_1 M _1 (r) for y < 1. Bottom: Curves showing the behaviour of y'/y 
a function of y for y > 1 (the functions y \, y' 2 , y' 3 and versus y are defined 
eq. (11)). 






Fig. 2. Curves of v = V (r) / (GAfao) 1 ^ 4 as a function of s(= r/r^) and //. The 
functions /T are defined in eq. (11). In all graphs, = 3/4, u = 1/4, and M = Mq, 
where Mo is an arbitrary mass. The functions ad and u are defined as ad = Md/M 
and u = r s /rd, where Mj is the mass of the disc and r s () is the effective radius 
of the spherical bulge (disc). 



Fig. 3. Curves of v = V (r) / (GAfao) 1 ^ 4 as a function of s (= r/r^) and ji'. The 
functions ji' are dehned in eq. (11). The solid line is the curve obtained using Mil- 
grorn’s original proposal, eq. (5). For all curves M = Mq, where Mq is an arbitrary 
mass, ad = 3/4, u = 1/4. The functions ad and u are defined as ad = Md/M and 
u = r s /rd , where Md is the mass of the disc and r s () is the effective radius of the 
spherical bulge (disc). 
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